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Abstract 

We obtain the variational upper bound for the ground- state energy of two- 
dimensional antiferromagnetic Heisenberg model on a square lattice at arbitrary 
value of the anisotropy parameter using the two-dimensional generalization of 
Jordan- Wigner transformation. Our result can be considered as an upper bound 
for the perturbation theory series about the Ising limit. 

At present time two dimensional quantum spin systems attract much attention in 
connection with the problem of high-Tc superconductivity. For the antiferromagnetic 
Heisenberg model at some values of the anisotropy parameter the existence of the long- 
range order was proved however the exact ground state is not known. Apart from 
the linear spin wave theory 0] various methods to evaluate the ground state energy for 
the Heisenberg antiferromagnet were proposed. For instance the perturbation theory 
and the cluster expansion about the Ising limit were used 0. However, although the 
convergence of the series of the perturbation theory is good these estimates are not 
the variational ones. At the same time the energy corresponding to any reasonable 
variational ground-state wave function cannot be computed exactly (for example of 
these calculations see ref. 0]). Finally at present time the accuracy of the numerical 
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simulations ||^ is not sufficiently high. In this context the variational estimates of the 
ground state energy for the two-dimensional antiferromagnetic Heisenberg model are 
of interest. 

In the present letter we obtain the exact upper bound for the ground-state energy of 
the s = 1/2 quantum antiferromagnet for arbitrary value of the anisotropy parameter. 
Our variational estimates are sufficiently low and may be useful in connection with the 
study of two-dimensional spin models in the framework of the other methods. 

Our method is based on the transformation which change the statistics of particles 
on a two- dimensional lattice. There are several ways to define the Hamiltonian of 
particles obeying the fractional statistics (anyons) on a lattice (for example see ref. [0). 
We can use the most natural form of the definition of anyon operators in terms of the 
fermions 

bf{a) = af exp -ia ^ 0i;n; , n; = a+a^ (1) 

where the operators af,ai obeys Fermi statistics and (pu is the angle between the 
direction from the site i to the site / and some fixed direction, x- axis for example. In 
accordance with the multi-valuedness of the anyon wave function the operator bf{a) is 
multi-valued at arbitrary fractional value of the statistical parameter a, which describes 
in Eq.(0) the deviation from the Fermi statistics. In particular, at a = 1 the operators 
(|l]) are the hard core boson operators, which commute at different sites and behave 
like the fermions at the same site. Expressing the spin operators {s = 1/2) in terms of 
the Holstein - Primakoff bose operators 

St = bt, Sr=b., S! = bth-^, 

we obtain the representation of spin operators in terms of the fermions which can be 
thought of as a two dimensional generalization of the well known Jordan - Wigner 
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transformation for one dimension: 

S+ = al exp (^-i ^ (f)un}j , S- = ajai - ^. 

The Hamihonian of the Heisenberg antiferromagnet H = J2<ij> SiSj has a comphcated 
form 

H = -^Yl I ^t^j exp I -i ^iji^i 1 + ^-c- ] + hi^j - (2) 

^ <ij> V V ly^td / J <ij> ^ ^ 

where (^iji = (pu — (pji and < ij > denotes the nearest neighbour sites. The minus sign 
before the first term in Eq.(^ is due to the redefinition of the operators Oj —Qi 
on one of the sublattices of the square lattice. In order to simphfy the Hamiltonian 
one can make the substitution rti fii in the exponential of Eq.@, where fii is the 
average particle number at a given site, n = 1/2 for the half filling [S^ = 0), which will 
be considered in the present paper. This procedure is usually refered to as a (vector) 
mean field (MF) approximation. After this substitution the MF Hamiltonian 

Hmf = -^Y1 iXijafaj + h.c.) + U ^ (n, - ^){nj - (3) 

<ij> <ij> 

where Xij = exp(— i X); ^ii/?^;)) describes the system of fermions in the homogeneous 
statistical magnetic field with the magnitude corresponding to the flux (p = it through 
the plaquette. The parameter U = 1 for the isotropic model. Due to the gauge 
invariance (a^ ajexp(^j)) the phases of Xij are depend on the gauge fixing condition. 
The sum of the phases around the closed contour is fixed and equal to the one-half 
flux quantum through the plaquette for the case of the half filling. The eigenstates of 
the MF Hamiltonian does not depend on the choice of the gauge. The second term in 
Eq.(|^) is the interaction of fermions. 

We use the variational theorem proved in ref.|0] for the hard core bosons in the 
absence of the interaction term. Let ipMpih, ■ ■ - iN) and Emf to be respectively the 



exact ground-state wave function and the ground state energy of the MF Hamiltonian 
(H) {ii, . . An - are the particle positions, < 4'mf\4'mf >= !)• Consider the contribution 
of a given bond to the expectation value of Eq.@ over the ground state. We have the 
following inequality: 

- \'^MF{i,i2, ■ ■ ■iN)\\i^MFUA2, ■ ■ ■iN)\ 

i2...ijv 

- ~^^[Y1 i^*MF{hi2,---iN)lpMF{j,i2,---iN)\- 

The left-hand side of this inequality is the contribution to the expectation value of the 
exact Hamiltonian in the bosonic representation. The normalization as well as the 
expectation value of the operator given by the last term of Eq. (|^) are the same for the 
wave functions "^mf and Thus it is proved that the ground- state energy of 

the initial bosonic Hamiltonian Eq is bounded from above by Emf'- 

Eq < Emf- 

This relation allows one to obtain an upper bound for the energy of the antiferromagnet. 
We have to obtain the appropriate variational estimate for the ground- state energy 
of the MF Hamiltonian (^. As a variational wave function let us choose the wave 
function corresponding to the Hamiltonian which is obtained from Hmf in the mean 
field approximation in respect to the fermion interaction. We assume the existence 
of Neel order in this state. Linearising the interaction and using the substitution 
< rii >^ (— l)*A/4 (we use the notation (—1)* = (—1)*^+*!') we obtain the Hamiltonian 

- ^ E + h.c.) - A Y.i-^yn^. (4) 

<ij> i 

In this formula A is the variational parameter which is to be determined from the 
condition of minimum of the expectation value of Hmf in the state given by the ground 
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state of the mean field Hamiltonian Eq.(^. This expectation value is the variational 
bound for the energy Emf- Note that the choice of the wave function is consistent with 
the MF treatment of the statistical interaction since the sum of the phases around the 
plaquette for the Neel ordered state is the same as in the case of = 1/2. 
The calculations are most easily performed using the symmetric gauge 

1 . 1 . 

Xi,i+x = ^(1 + ^(-1)'), Xi,i+y = ^(1 - ^(-1)'), 

where x, y are the unit vectors corresponding to the lattice spacing. In the momentum 
space in terms of the doublets ipik = {o,k,cik^Q) (0 < kx,ky < it, Q = (vr,7r)) and 
i'2k = {ak 

a^-gj (0 < —kx, ky < TT, Qi = (— vr, tt)) the equation (^) has the form 
H= J2 ^tkMk^Pik+ E i^tkMk^p2k, 

kx>0, ky>0 

where the matrix is 

/ 



Ci A — ?C2 I 1 

Cl 2 = —F^ (cos kj; ± COS ky). 

V2 



\ A + ic2 -Cl 

The eigenvalues are Ek = ±(cos^ A;^ + cos^ ky + A^)^/^ where the momentum k is 
restricted to the half of the Broullien zone ky > 0. The negative energy levels are 
filled. Let us calculate the average of Hmf over this state. The values of < Xij^t'^j > 
for a given < ij > does not depend on the choice of the gauge. This values are real 
(and positive) which can be deduced from the parity invariance of our state. The 
expectation value of the second term in Eq.(^) is < n^rij >=< rii >< rij > — < 
afttj >< a^ai >. The expression for the particle number at a given site has the 
form < rii >= l/2 + (— l)*Ai/4, where the parameter Ai does not coincide with the 
parameter A. We obtain for Ai and ^ = | < afaj > \ (which is the same for all bonds) 
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the following expressions: 



A 



E 



COS' 



i kx ~\~ COS k', 



'y 



kx>0, ky>0 



kx>0, ky>0 



Ek 



The final expression for the variational estimate is 



16f/ 



A2 



fca;>0, ky>Q 



E 
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(5) 



where is the number of sites. The sum of the first two terms is the energy in the 
mean field approximation Eq.(^. 

At f/ = we get the exact energy of the MF Hamiltonian and the corresponding 
estimate for the energy for the XY model is —0.2395 per bond. In comparison with 
the energy determined with the help of the numerical simulations — 0.27(±10%) 0, 
the bound is too high. It is less restrictive than the bound based on the simple trial 
variational wave function wave function. For example the energy corresponding to 
the Neel ordered state (in the y-direction) is —0.25. That is in agreement with the 
statement of ref's.[0, ^ that the corrections due to the fiuctuations around the average 
magnetic field background are of order of unity. The situation is different for the 
isotropic (XXX) model {U =1). In this case the perturbation theory series is rapidly 
converges and the corrections due to the statistics of particles are suppressed. For 
instance for the Hamiltonian (0) the corrections to the MF approximation are of order 
~ l/(2f/)^ which is a sufficiently small value |]^. In this sense our result can be 
considered as an estimate from above for the perturbation theory series about the 
Ising limit. It is difficult to establish the restrictions of this type using the other 
methods. Minimizing the expression (^ with respect to A (Aq = 1.19) we obtain 
£^^^^/2L^ = —0.33034, which is sufficiently good upper bound for the energy. For 
comparison the best estimate obtained using the method of ref.0 is —0.334. Note 
that although the prediction of the linear spin wave theory |^ is —0.329, this method 
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does not result in the correct ground-state wave function and this value cannot be 
considered variational bound. 

For the anisotropic model we proceed as follows. For simplicity let us consider the 
axially symmetric model although our method can be easily generalized to the case 
of arbitrary asymmetry. We use the description in terms of the Holstein-Primakoff 
bosons for the equivalent Hamiltonian H = X]<«j> 

(5f 5J + ^SfS] + After the 

substitution bi (— we get 

H = -Y: C-^iKbj + h.c) + ^(6+6+ + h.c)) + E (^^ - ^)(^. - \)- (6) 
<ij> ^ ^ ^ ^ <ij> ^ ^ 

Consider the trial variational wave function with the fixed number of bosons. For 
this state the expectation value of the second term ~ ibitb^ + bihj) in Eq.(|^) is zero. 
The Hamiltonian without this term can be used to obtain the variational estimate 
for the ground-state energy for arbitrary 7 according to our method. Note that the 
contribution of the omitted term is small in the framework of the perturbation theory 
[0] since it appears only in the fourth order. The analysis can be performed at arbitrary 
value of the parameter 7. For the XY model (7 = 0) which is equivalent to the system 
of the hard core bosons at the half filling we found the estimate —0.26776 per bond 
(Ao = 3.4). This estimate is in agreement with result of the numerical simulations 0. 

In conclusion, although the wave function corresponding to the mean field Hamilto- 
nian cannot be used to describe the long-range properties of the model (for example, 
the energy of the low-lying excitations) the ground-state energy can be estimated with 
the sufficiently high accuracy. We found the variational upper bound for the ground- 
state energy of the two-dimensional Heisenberg antiferromagnet on a square lattice 
at arbitrary value of the anisotropy parameter. Our results can be thought of as a 
peculiar upper bound for the perturbation theory series about the Ising limit and may 
be useful in the context of the other approaches. 
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